FORMPTO.1390 
(REV 11-2000) 



JC05R8e'dPCT/PTO 1 4 MAR 2002 



TRANSMITTAL LETTER TO THE UNITED STATES 

DESIGNATED/ELECTED OFFICE (DO/EO/US) 
CONCERNING A FILING UNDER 35 U.S.C. 371 



ATTORNEY DOCKET NUMBER 

0887-41 47US2 



U.S. APPUCATIOSJ NO. (If known, see 37 CFR 1.51 

IO/Ot8a8190 



INTERNATIONAL APPLICATION NO. 

PCT/USOO/25136 



INTERNATIONAL FILING DATE 
14 September 00 (14.09.00) 



EARLIEST PRIORITY DATE CLAIMED 
14 September 99 (14.09.99) 



TITLE OF INVENTION 
METHOD AND SYSTEM FOR IMAGING THE DYNAMICS OF SCATTERING MEDIUM 



APPLICANT(S) FOR DO/EO/US 
Randall L. BARBOUR 



Applicant herewith submits to the United States Designatecl^lected Office (DO/EO/US) the following items and other information: 

1. ^ This is a FIRST submission of items concerning a filing under 35 U.S.C. 371. 

2. □ This is a SECOND or SUBSEQUENT submission of items concerning a filing under 35 U.S.C. 371. 

3. ^ This is an express request to begin national examination procedures (35 U.S.C. 371(f)). The submission must include 

items (5), (6), (9) and 21 indicated below. 



6. 



7. 



^ The US has been elected by the expiration of 1 9 months from the priority date (Article 3 1 ). 
^ A copy of the International Application as filed (35 U.S.C. 371(c)(2)) 

a. □ is attached hereto (required only if not communicated by the International Bureau). 

b. g| has been communicated by the Intemational Bureau. 

c. □ is not required, as the application was filed in the United States Receiving Office (RO/US). 

□ An English language translation of the Intemational application as filed (35 U.S.C. 371(c)(2)). 

a. Q is attached hereto. 

b. □ has been previously submitted under 35 U.S.C. 154(d)(4). 

^ Amendments to the claims of the Intemational Application under PCX Article 19 (35 U.S.C. 371(c)(3)) 

a. □ are attached hereto (required only if not communicated by the Intemational Bureau). 

b. □ have been communicated by the hitemational Bureau. 

c. □ have not been made; however, the time hmit for making such amendments has NOT exphed. 

d. 13 have not been made and will not be made. 

□ An English language translation of the amendments to the claims imder PCT Article 19 (35 U.S.C. 371(c)(3)). 
[3 An oath or declaration of the inventoi(s) (35 U.S.C. 371(c)(4)). 



10. □ An EngUsh language translation of the annexes to the Intemational Prelimmary Examination Report under PCT 
Article 36 (35 U.S.C. 371(c)(5)). 

Items 11. to 16. below concern document(s) or information included: 

1 1. El An Information Disclosure Statement under 37 CFR 1.97 and 1.98 with Form PTO 1498 and copy of cited reference. 

12. □ An assignment document for recording. A separate cover sheet m compliance with 37 CFR 3.28 and 3.3 1 is included. 

13. □ A FIRST preliminary amendment. 

14 □ A SECOND or SUBSEQUENT preliminary amendment. 

15. □ A substitute specification. 

16. □ A change of power of attomey and/or address letter. 

17. □ A computer-readable form of the sequence listing in accordance with PCT Rule 13ter.2 and 35 U.S.C. 154(d)(2) with 
copy of Statement Under 37 CFR Section 1.821(f) and WIPO Standard ST.25 as filed with the Intemational Bureau 
of WIPO. 

18. □ A second copy of the published intemational application under 35 U.S.C. 154(d)(4). 
□ ^ second copy of the English language translation of the intemational application under 35 U.S.C. 154(d)(4). 



687209 vl 



^ISnec'dPCT/PTO 1 4 MAR2QQ2 



APPLICATION NO. (if known, see 37 CF.R. 1.51) 

10/088190 



INTERNATIONAL APPLICATION NO. 

PCT/USOO/25136 



ATTORNEY'S DOCKET NO. 

0887-4147US2 



20. S Other items or information: 

- Copy of published PCX Application - Publication number WO 01/20305, published on 22 March 01 with International Search Report. 



m 
ra 



687209 vl 



^BmdPmim t 4 MAR 2002 



U.S. APPLICATION NO (if known, see 37 C F R. I 51 

10/ft88190 



INTERNATIONAL APPLICATION NO 

. PCT/USOO/25136 



ATTORNEY'S DOCKET NUMBER 

0887-4147US2 



2L K| The following fees are submitted: 




CALCULATIONS PTO USE ONLY 


BASIC NATIONAL FEE (37 CFR 1 .492 (a) (1) - (5) ): 

Neither international preliminaiy examination fee (37 CFR 1.482) 

nor international search fee (37 CFR 1.445(a)(2) paid to USPTO 

and International Search Report not prepared by the EPO or JPO $ 1 040.00 




International preliminary examination fee (37 CFR 1.482) not paid to 

USPTO but International Search Report prepared by the EPO or JPO $890.00 




Tntf^matimnal nrplimttiarv eviimififltion flsft C^l {^PT? 1 48'!?^ nnt nfliH tn 

USPTO but international search fee (37 CFR 1 .445(a)(2) paid to USPTO $740.00 




International preliminary examination fee paid to USPTO (37 CFR L482) 

but all claims did not satisfv orovisions of PCT Article y^iW- $690 00 




International preliminaiy examination fee paid to USPTO (37 CFR 1.482) 

and all claims satisfied provisions of PCT Article 33(1) - (4) $100.00 




ENTER APPROPRIATE BASIC FEE AMOUNT = 


$690.00 




Surcharge of$130 for furnishing the oath or declaration later than □ 20 □ 30 
months from the earliest claimed priority date (37 CFR L492(e)). 






CLAIMS 


NUMBER FILED 


NUMBER EXTRA 


RATE 




Total claims 


25-20 = 


5 


X $18.00 


$90.00 




Independent claims 


2-3 = 


0 


X $84.00 


$0.00 




MULTIPLE DEPENDENT CLAIM(S) (if applicable) 


$280.00 


$0.00 




TOTAL OF ABOVE CALCULATIONS = 


$780.00 




^ Applicant claims small entity status. See 37 CFR. 1 .27. The fees indcated above 
are reduced by 1/2. 


$390.00 




SUBTOTAL = 


$390.00 




Processing fee of $130.00 for furnishing the English translation later than □ 20 □ 30 
months finm the earliest claimed priority date (37 CFR 1.492(f)). 


$ 




TOTAL NATIONAL FEE = 


$390.00 




Fee for recording the enclosed assignment (37 CFR 1.21(h)). The assigranent must be 
accompanied by an ^propriate cover sheet (37 CFR 3.28, 3.31). $40.00 per property X — 


$ 




TOTAL FEES ENCLOSED = 


$390.00 












Amount to be 
refunded: 


$ 










charged: 


$ 



a. 13 A check in the amount of $390.00 to cover the above fees is enclosed. 

b. □ Please charge my Deposit Account No. 13-4500, ORDER NO. m the amount of $0. 00 to cover the above fees. 

c. E] The Commissioner is hereby authorized to charge any additional fees which may be required, or credit any 

overpayment to Deposit Account No. 13-4500, ORDER NO. 0887-4 147US2. A duplicate copy of this sheet is enclosed. 

d. □ Fees are to be charged to a credit card. WARNING: Information on this form may become public. Credit card 

information should not be included on this form. Provide credit card information and authorization on PTO-2038. 
NOTE: Where an appropriate time limit under 37 CFR 1.494 or 1^195 has not been met, a petition to revive (37 CFR 
1.137(a) or (b)) must be fUed and granted to restore the application to pending status. 

SEND ALL CORRESPONDENCE TO: 

Morgan & Finnegan LLP 
345 Park Avenue 
New York, NY 10154-0063 
Telephone: 212-758-4800 
Telecopier. 21 2-751 -6849 




Registration No. 41,190 



687209 v1 



10/088190 

\ %l JC13 Rec'd PCT/PTO l 4 MAR 2002 



wo 01/20305 IP^'*^ PCT/US00/2S136 

METHOD AND SYSTEM FOR IMAGING THE DYNAMICS 
OF SCATTERING MEDIUM 

This invention was made with U.S. Govenunent support under contract number 

5 CA-RO166184-02A, awarded by the National Cancer Institute. The U.S. Govenunent 

has cotain lights in the invention. 

This ^plication claims tiie benefit under 35 U.S.C. §120 of prior U.S. Provisional 

Patent Application Serial Nos. 60/153,926 filed Septembo- 14, 1999, entitled DYNAMIC 

TOMOGRAPHY IN A SCATTERING MEDIUM and 60/154,099 filed September 15, 

f-f 10 1 999, entitled DYNAMIC TOMOGRAPHY IN A SCATTERING MEDIUM. 

111 This application is related to copending application serial number "not vet 

i ■ 

C3 assigned", attorney docket number 0887-4147PC1, filed on the same date as this 

application, entitled "SYSTEM AND METHOD FOR TOMOGRAPHIC IMAGING OF 
DYNAMIC PROPERTIES OF A SCATTERING MEDIUM" by inventors R. Barbour 
W 15 and C. Schmitz and is hereby incorporated by reference (hereinafter the **Barbour 
4147PC1 application**). 

This application is also related to copending application serial number "not vet 
assigned" attorney docket number 0887-41 49PC2, filed on the same date as this 
application, entitled "IMAGING OF SCATTERING MEDIA USING RELATIVE 
20 DETECTOR VALUES" by inventor R. Barbour and is hereby incorporated by reference 
(hereinafter the "Barbour 4149PC2 application")^ 

Field of the Invention 

The invention relates to the field of imaging in a scattering medium, and in 
particular to methods of imaging the dynamics of a scattering medium. 

1 
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Background 

Imaging in a scattering medium relates generally to the methods and techniques 
for generating an image of the internal properties of a scattering medium based on the 
detection of scattered energy. 

Many systems and techniques have been developed for imaging of scattering 
media. A typical system for imaging based on scattered energy detection includes a 
source for directing energy into a target medium and a plurality of detectors for 
measuring the scattered energy exiting the target medium at various locations with 
respect to the source. Based on the measured energy exitmg the target medium, it is 
possible to reconstruct an image representation of the cross*sectionaI scattering, 
absorption or other properties of the target medium. The values of the absorption and 
scattering properties of the medium can vary depending on the wavelength and types of 
energy employed as an imaging source. These values are also frequently spatially 
varying. These techniques permit the use of types of energy and wavelengths, such as 
near infrared light energy, that are not suitable for projection imaging techniques, such as 
x-ray imaging. Thus these techniques have great potential for detecting and imaging 
properties of media, such as human tissue, that can not be revealed using energy sources 
commonly employed in projection imaging methods. 

Exemplary methods and systems for imaging of scattering media are disclosed in 
Barbour et al, U.S. Patent No. 5,137,355, entitled "Method of Imaging a Random 
Medium," (heremafter the "Barbour *355 patent"), Barbour, US. Patent No. 6,081,322, 
entitled **NIR Clinical Opti-Scan System," (hereinafter flie "Baibour *322 patent"), the 
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Barbour 4147PC1 j^jplication, the Barbour 4149PC1 application and the Baibour 
4147PC2 application. 

As can readily be appreciated, there are many instances where use of these 
techniques are highly desirable. For example, one flourishing application area is in the 
field of optical tomography. Optical tomography typically uses near infrared (NIK) 
energy as an unaging source. Contrary to imaging methods relying on the use of ionizing 
radiation and/or toxic/radioactive contrast agents, NIR optical tomographic imaging 
methods bear no risk of causing harm to the patient. The dose of optical intensity used 
remains far below the threshold of thermal damage and is therefore safe. In the regime of 
wavelength / intensity / power used, there are no effects on the tissue that accumulate 
with increasing light energy dose due to over-all irradiation time. 

Other favorable attributes of optical tomography include the use of low*cost, 
potentially portable devices that employ highly integrated, economical off-the-shelf data 
processing electronics and semiconductor lasers {laser diodes). Such features contrast 
with other imaging technologies commonly used in clinical diagnosis that require large, 
fixed facilities such as MRI and x-ray CT imaging. Additionally, since a significant 
computational effort may be required for both image reconstruction and data analysis, the 
technology particularly gains from the exponential growth in the ratio of computmg 
power to cost. 

It is well appreciated that optical tomography has the potential to provide insights 
into anatomy and physiology that are unavailable from other imaging methods. For 
example, optical tomography, using near infrared energy, can identify the spatial 
variations in blood volume and blood oxygenation levels because of its sensitivity to 
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hemoglobin states. These measures have considerable potential value in diagnosing a 
broad range of disease processes that are known to influrace hemoglobin states. 

For example, a common feature of breast tumors, and solid tumors generally, is 
the occurrence of neovascularized tissue. Ultrastructurally, these tissues are highly 
5 disorganized and exhibit functional abnormalities. Often the microvessels are dilated, 
tortuous, elongated and saccular. There is excessive branching of the vessels, including 
significant arterio-venous shunting as well as blind vascular endings. Aberrant vascular 
morphology and decreased vessel density are responsible for mcrease resistance to flow. 
Q The resistance to flow combined with an enlarged diffusion distance, due to the 

10 expansion of the extravascular space, can lead to perfusion with hypoxemic and nutrient- 
deprived blood. The net effect of this state is the occurrence of substantial spatial and 
temporal heterogeneity in the tumor metabolic microenvironmrait, 

Ijfl Although these attributes of disease tissue are well appreciated, the availability of 

■ a suitable detection methodology able to take full advantage of these characteristics has 

fS 

p j 1 5 been notably lacking. An appropriate methodology would be one sensitive not only to 

altered hemoglobin states (i.e., localized variations in tissue blood volume and 

oxygenation states), but also to their dynamics under homeostatic conditions or in 

response to specific provocations. 

A variety of methods involving imaging and non-imaging modalities are available 

20 for assessing specific features of the vasculature. Detailed images of the vascular 

architecture involving larger vessels (> 1 mm dia.) can be provided using x-ray enhanced 

contrast imaging or MR angiography. These methods however are insensitive to 

hemoglobin states and only indirectly provide measures of altered blood flow. The latter 

4 
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is well accomplished, in the case of larger vessels, using Doppler ultrasound, and for 
near-surface microvessels by laser Doppler measurements, but each is insensitive to 
variations in tissue blood volume or blood oxygenation. Ultrasound measurements are 
also limited by their inability to penetrate bone. 

In principle, imaging methods based on the detection and analysis of scattered 
energy, such as optical tomograpic methods, can provide either direct or indirect 
measures of all of these parameters. However, the known methods and systems have 
several shortcomings. First, known methods and systems provide images having low 
contrast and resolution. Second, these me&ods and systems do not image the dynamic 
properties of highly scattering media. Third, these methods and systems require accurate 
calibration and are susceptible to errors. There are several reasons for these problems 
with known systems and methods. These reasons relate principally to how measurements 
are performed and how measurement data is analyzed. 

For example, when imaging human tissue, the natural occurrence of vascular 
frequencies arising from cardiac, respiratory and vasomotor activity, produces time 
variations in, for example, the absorption properties of tissue due to changes in tissue 
blood volume. Significantly, the process of vasomotion, perfusion first in one region, 
then another, can be expected to produce spatially convolved images should data be 
collected on a time frame that is long compared to the reciprocal of the frequency of these 
processes. Thus methods that collect time-averaged data will predictably yield images 
whose contrast and resolution are degraded by such variability. 

Also influencing the quality of reconstructed images, is the approach used to 
analyze the acquired data. Many of the known imaging schemes consider, in some 
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manner, the comparison of measured values to predicted values. Typically, these 
methods, mcluding that described in the Barbour '355 patent, employ numerical methods 
that seek to minimize the difference between sets of measured and predicted values, and 
in doing so seek to provide improved estimates of the properties of the unknown target 
medium. These analysis schemes, referred to as model based methods, assume 
equivalency in the efficiency of measured detector values and computed predicted values. 
Although the derivation of accurate estimates of the efficiency of measured responses is 
possible in principle, in practice, the natural plasticity of tissue, its mainly arbitrary shape 
and variable composition and noted variability in hemoglobin levels, all serve to 
confound efforts to devise practical methods that provide reliable estimates. 

In addition, the physics of light transport in highly scattering media, such as 
tissue, imposes further practical constraints that relate to the method adopted for data 
analysis. At issue is the required accuracy of assumptions made in order to generate 
predicted detector values, especially those adopted for the initial estimate. These 
assumptions are commonly referred to as the "initial guess". Small errors in the initial 
guess of optical properties of the reference medium can lead to large errors in the 
computed detector values. One consequence of this can be the severe corruption of the 
information content of the data vectors leading to artifact-laden images. Adding to the 
mentioned uncertainties is the well-known property of reconstruction methods that 
employ Unear operators regarding their sensitivity to undetermined data sets {i,e., 
insufficient amount of collected data) and measurements based on restricted views (e.g., 
backscatter only, transmission only). The net effect of these limitations can render 
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solutions to image recovery problems of this type overly sensitive to the influence of 
experimental noise (i.e.» ill-conditioned), provide nonunique solutions (ill-posed) or both. 

These concerns are well appreciated by those skilled in the art of image 
reconstruction methods. It is also well appreciated that, in general, there are no simple or 
S well-defined methods diat can be universally applied to overcome the noted limitations. 
In this regard, specification of suitable conditions that can satisfactorily deal with the 
noted concerns is an art whose successful implementation requires considerable skill. 

In addition to the need to provide stable solutions to the image reconstruction 
problem^ consideration of the infonnation content of the reconstructed image has 
10 considerable importance. As presently practiced, the method of optical tomogr^hy 
considers the evaluation of static states or employs time averaging methods to minimize 
the influence of signal instability originating fi-om tissue dynamics. 

The goal of these studies is to provide image maps that define spatial variations in 
the optical properties of tissue (usually, absorption and scatter) from which may be 
IS d^ved, for instance, estimates of spatial varying hemoglobin states. It is understood 
however, that the latter is fundamentally governed by dynamic processes whose details 
have the potential to reveal a wealth of infonnation regarding functional features of the 
vasculature, in particular as it relates to its interaction with surrounding tissue. 

The ability to measure dynamic processes of a medium can reveal information 
20 that is unobservable from static or time-averaged measures. In the case of physiological 
systems, the form of the dynamic process has added significance. For instance, it is well 
understood that many time varying processing have an underlying nonlinear character. 
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Nonlinear dynamic processes, in biological systems, are often chaotic and exhibit the 
characteristic feature of sensitivity to initial conditions. 

The existence of such behavior has important implicatioi^ in the understanding of 
disease processes and well as for (he ^proacbes taken for therapy. For instance, the 
approach needed to control a chaotic system is quite different from that for a linear 
system, wherein the system response is proportional to the magnitude of the input 
stimulus. Thus it has been proposed that more effective therapies can be realized from a 
series of well-timed perturbations rather than from the standard approach of applying a 
constant stimulus, the method commonly used in many pharmacological interventions. 
Also of interest, and related to this, is the seemingly general finding that the occurrence 
of chaotic behavior in physiological systems is a sign of health and its absence is a sign 
of disease. 

For instance, it is known that heart rate variabihty is chaotic. Significantly, loss 
of this signature with the appearance of periodic oscillations is among the strongest 
predictors of sudden cardiac death. A similar phenomenology has been observed in 
infants who succumb to sudden infant death syndrome. In this case, the normally chaotic 
respiratoiy rate becomes periodic prior to the fatal incident. Similarly, during epileptic 
seizures, electroencephalographic recordings exhibit a transition from chaotic to periodic 
activity. 

Presently, the capacity to monitor dynamic behavior in vascular structures is 
limited principally to near surface measures using laser Doppler methods. Measures of 
the time variability of the vascular caliber and flow motion for larger vessels is possible 
using Duplex ultrasound. However, these measures are insensitive to the activity of the 

8 
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microvasculature, do not provide for foil cross-sectional views, and are not sensitive to 
the dynamics of hemoglobin states. 

Although optical methods, such as Laser Doppler, pulse oximetry, 
photoplethysmography and the like, can be used to monitor dynamic states of the 
vasulature, none are capable of providing such measures in the form of a cross-sectional 
image, especially in the case of large tissue structures. Moreover, known optical methods 
for cross-sectional imaging of the properties of a scattering media have not bem used to 
derive dynamic measures of these states, and are plagued with a host of technical 
limitations, such as low contrast, low resolution images, prolonged computing times, 
excessive sensitivity to errors in initial estimates. These limitations render it unlikely that 
usefol information regarding dynamic states could be derived from these known methods. 

Overcoming the indicated drawbacks and concerns is critical for widespread 
practical implementation of optical tomography as a diagnostic tool because (1) improved 
contrast and resolution are essential to feature identification and visualization, (2) a static 
(snapshot) or time averaged image of a time evolving property does not provide 
discovery of the physiological dynanwc processes and (3) measures of dynamic processes 
can yield critical information needed for improved diagnostic methods and therapies. 

For the foregoing reasons, there is a need for a method of improvmg the contrast 
and resolution of reconstructed images. There is also a need for a method of imaging 
dynamic properties of dense scattering media, especially as it relates to dynamic 
prop^ies of vascular states in large tissue structures as revealed by time variation in 
hemoglobin states. There is yet a further need for a method that can provide dynamic 
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images without undue reliance on complex calibration schemes or computationally 
intensive ntmierical methods. 

SUMMARY OF THE INVENTION 

The present invention satisfies these needs by providing a method for (1) 
enhancing resolution and contrast of various dynamic features of the medium being 
imaged, (2) imaging the dynamic properties of the target medium, and (3) generating 
images of the dynamic properties using techniques that reduce the need for undue system 
calibration and produce more stable solutions with the reconstruction algorithm. 

It is one object of the invention to generate a map of the dynamic properties of a 
scattering medium in a cross-sectional view from a time series of collected data 
measurements. The measurements are obtained by directing energy into the target during 
a period of time and measuring the energy emerging from the target during the period of 
time, whoreby a time series of measured energy is collected. 

It is another aspect of the invention to generate a map of the cross-sectional 
dynamic properties of the target medium in a cross-sectional view from a time series of 
images of the cross-sectional propoties of the target medium, wherein the dynamic 
properties are extracted from the time series of images using time series analysis 
methods. 

It is a further aspect of the invention to generate a map of the cross-sectional 
dynamic properties of the target medium from the time series of measurements at each 
detector, wherein the dynamic properties are extracted from the time series of 
measurements at each detector using time series analysis methods. 
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It is a further aspect of the invention to construct a map of the dynamic properties 
of the medium, wherein the measured energy is processed using a modified perturbation 
fonnulation referred to as the normalized difference method, based on the radiation 
transport equation, using relative energy measurements. The relative energy 
measurements being, for example, the relative difference between an instantaneous 
measure and a time average mean of measures. 

BRIEF DESCRIPTION OF THE HGURES 

For a better understanding of the invention, together with the various features and 
advantages thereof reference should be made to the following detailed description of the 
preferred embodiments and to the accompanying drawings wherein: 

FIG. 1 is a schematic illustration of an exemplary imaging system; 

FIG. 2A is an illustration of a mechanism used to rhythmically inflate two 
balloons; 

FIG. 2B is a plan view of the illustration of FIG. 2A; 

FIG. 2C is an illustration of the geometric arrangement of optical fibers; 

FIG. 3 is a reconstructed image of the cross-sectional properties of a target 
medium shown in FIG. 2A at an instant in time; 

FIG. 4 is a cross-sectional map of the 0.2 Hz component of the amplitude of the 
image sequence's DFT for the target medium shown in FIG. 2; 

FIG. 5 is a cross-sectional map of the 0.2 Hz component of the phase of the image 
sequence's DFT for the target medium shown in FIG. 2 with no phase difference between 
the oscillating balloons; 

11 
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FIG. 6 is a cross-sectional map of the 0.2 Hz component of the phase of the image 
sequence's DFT, 180^ phase difference between oscillating balloons; 

FIG. 7 A is a schematic illustration of apparatus with two-balloons; 

FIG. 7B is a series of reconstructed cross-sectional images Ifrom DFT analysis of 
image time series at 0.1 Hz and 0.24 Hz; 

FIG. 8 is a cross-sectional map of ratio of FT amplitude of cardiac to respiratory 
frequency from foreann study; 

FIG. 9 is an reconstructed MR image of foreann, with identified anatomical 
structures; 

FIG. 10 is an overlay of images shown in FIGS. 8 and 9; 

FIG. 1 1 A is a graph of the CSD spectrum amplitude at position (12,19); 

FIG. I IB is a graph of the CSD spectrum amplitude at position (1 1,10); 

FIG. 1 1 C is a graph of the CSD spectrum amplitude at position (1 8,27); 

FIG. 12 is a cross-sectional map of amplitude, xlO*, of the DFT at the finger-flex 
frequency; 

FIG. 13 is an overlay of the map shown in FIG. 12 and an anatomical image 
(MRI) of the same tissue structures; 

FIG. 14 is a graph of the temporal variations in //^ (absoiption coefficient) for 
pixels located within the involved muscle groups; 
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FIG. 15 is a graph of the time-dependent intensity variations for foreann 
measurements observed for a detector positions opposite the source; 

FIG. 1 6 is a cross-sectional map of the correlation dimension, v, computed, in 
each pixel, from time series of reconstructed images of (810 nm) obtained from 
5 representative detector data shown in FIG. 1 5; 

FIG. 17 is a cross-sectional map of the correlation dimension, v, computed in each 
pixel from time series of reconstructed images of /j^ (8 1 0 nm) obtained from 
measurements performed on a solid white Dehin® rod; 



|j FIG. 18 is a cross-sectional map of (largest value of Lyapunov exponent) in 

^ ..' , 10 each pixel derived from the same time series of reconstructed /4 (810 nm) images as used 



to generate v map shown in FIG. 1 6; and 



4^ 

£3 FIG. 19 is a cross-sectional comparison map of the m map in FIG. 18 to 

M corresponding results obtained from surrogate data sets. 



DETAILED DESCRIPTION OF THE INVENTION 

IS Introduction 

The present invention is built upon known basic approaches to tomographic 
imaging of a scattering medium, by enabling the investigation of the spatio temporal 
dynamics of a target medium. This enhancement represents a fundamentally new 
imaging modality with broad based applications such as in industrial processes involving 
20 preparation of powdered mixtures, imaging subsurface turbulence in murky waters and 
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throughout clinical medicine both as a highly sensitive diagnostic imaging tool, as well as 
for use in treatment planning, therapeutic monitoring and follow-up. 

For example, the present invention recognizes that the occurrence of dynamic 
behavior in human tissue lends the collected data and reconstructed images to evaluation 
using time series analysis methods, including pattern recognition schemes. Dynamic 
behavior is known to occur in vascular structures in response to cardiac, respiratory, 
vasomotor and local metabolic influences. Because optical tomography, based on near 
infrared measurements is sensitive to hemoglobin, which is nearly always restricted to the 
vascular compartment time variations in the measured signal, can be used to reveal the 
vascular frequencies. Although, similar information can be obtained by 
photoplethysmography, it is not available for large tissue structures nor in an cross- 
sectional imaging modality* 

Detection of these temporal signatures/dynamic behaviors in a cross-sectional 
view provides a wealth of new information. For instance, the presence (or absence) of a 
cardiac frequency with respect to known anatomic landmarks (e.g., major arteries) 
provides a measure of the patency of localized arterial flow. Low-amplitude signals 
might suggest the presence of arterial stenosis. In addition, it is well appreciated that 
blood flow is under autonomic control. The absence or attenuation in amplitude of the 
respiratory or vasomotor frequencies might suggest an underlying peripheral neuropathy. 

Further insight concerning the presence or absence of basic control mechanisms 
can be obtained from time-frequency analysis methods. For instance, the low-frequency 
vasomotor response in the major arteries of the forearm are nearly synchronous, both 
ipsilaterally and contralaterally, and under tight autonomic control. Disorganization of 

14 
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this control might suggest the presence of pathology. In fact, loss of control mechanisms 
is a well known hallmark of neoplastic tissues. The nearly total absence of the vasomotor 
response in tumor tissue is of particular interest. 

Still another factor motivating dynamic measures is an appreciation of the 
5 information content available from the study of homeostatic challenges. It is well 
appreciated that the altered vascular architecture in tumors leads to a condition of 
sluggish perfusion. The existence of such states maybe be clearly revealed by induction 
of a homeostatic challenge. This could take on several forms. For example, simple 
manipulations such as deep breathing, response to cold, or mild venous congestion 
10 imposed by inflation of a pressure cuff, all have been shown to produce large amplitude 
local variations in the optical signals attributable to tissue blood volume in tiie forearm. 
Extension of these measures to the breast and other tissue could reveal the presence of 
p., tumors with high specificity by virtue of an altered hemodynamic response. 

There are other important implications to be considered concerning the value of 

f'^ 

2;! 1 5 measuring the time- varying properties of hemoglobin states in a cross-sectional view, 
' I For instance, the details of the measured dynamics at the various vascular frequencies 

need not be the same for measures of blood volume and blood oxygenation. In addition, 
given the known heterogeneity in tissue perfusion, even in healthy tissue, whatever 
dynamic features do exist will themselves be spatially var>dng. Moreover, altered 
20 dynamics in blood volume or blood oxygenation indicate altered local metabolic states, 
the existence of more central control deficits (e.g., autonomic, cardiac, respiratory), or 
both. Hence these measures can provide an assessment of an integrated physiological 
state, and also possess important features pertaining to tissue/vascular coupling. In fact, 
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all of the spatiotemporal features attributable to hemoglobin states that do exist will likely 
respond to a host of pharmacological agents and other treatment modalities. Such 
measures could serve to identify desired or tmdesired responses to therapy. 

Although these examples relate to the measure of hemoglobin states, it will be 
5 appreciated by one skilled in the art that there are a plethora of dynamic features that may 
be observable using tomographic measures in scattering media in human tissue as well as 
other scattering media, such as murky water, foggy atmospheres and the like. For 
instance nerve activation and muscle movement can also by measured owing to changes 
Q in scattering properties of tissue. In addition, various exogenous agents that produce 

10 contrast in absorption, scatter or fluorescence can be introduced to observe dynamic 



w4 

states of tissue. 

g| The present invention includes new data collection schemes and image analysis 

g| methods from which are extracted maps that reveal spatial temporal signatures of the 

target medium. These methods can include, linear and non-linear time series analysis 
1 5 methods as well as techniques for pattern recognition. 



"m 



There are three principal elements to practical dynamic imaging. The first 
element is the use of a fast, parallel, multi-channel acquisition system. For example, a 
system for use in dynamic optical tomographic imaging of tissue is described briefly 
below and in further detail in the copending Barbour 4147PC1 application. The second 
20 element is to evaluate the acquired tomographic data to produce a cross-sectional image 
that reveals perturbations of the reconstructed optical coefficients from, for example, a 
time-averaged mean. This perturbation method is described briefly below and in further 
detail in the Barbour 4149PC2 application. The third element is to collect a time series of 
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data and subject either the time series of measured data or a time series of reconstructed 
images from the measure data to analysis. This can include use of various linear and 
nonlinear time-series analysis methods, and related techniques such as pattern 
recognition, to extract features identifying dynamic behavior in the properties of the 
target medium. These methods are discussed in faribsr detail below. 

The System 

Numerous imaging systems such as those disclosed in the Barbour '355 patent, 
and the Barbour '322 patent have been developed for use in imaging of a scattering 
medium which are hereby incorporated by reference. 

A system providing high speed data capture of one or more wavelengths 
simultaneous is disclosed in the Barbour 4147PC1 application. This system is capable of 
capturing multiple wavelength data at rates up to 150 Hz and enables the reconstruction 
of cross-sectional images of real-time events associated with vascular reactivity in a 
variety of tissue structures (e.g., limbs, breast, head and neck). Fast data collection 
methods are particularly useful because there are many disease states with specific 
influences on the spatial-dynamic properties of vascular responses in hemoglobin states. 

A schematic illustration of an exemplary optical system for fast data collection is 
shown in FIG. 1. This system includes a computer 102, energy soxirces 104, 106, a 
source demultiplexer 108, an imaging head 110, detectors 1 12 and a data acquisition 
board 1 14. 

A target 1 16 placed in the imaging head 1 10 is exposed to optical energy from 

sources 104, 106. The optical energy originating from energy sources 104, 106, is 

combined by beam splitter 118 and is delivered to source demultiplexer 108. The source 
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demultiplexer 108 is controlled by computer 102 to direct the optical energy to source 
fibers 120 sequentially. Although two energy sources 1 04, 106 are shown in this 
embodiment, an additional number of energy sources, each having different wavelengths 
can be employed. Moreover a single variable wavelength energy source can be 
implemented such as Ti-Sapphire laser, dye lasers and the like. 

Each source fiber 120 carries the optical energy firom the demultiplexer 108 to the 
imaging head 110 where the optical energy is directed into (he target 1 16. The imaging 
head 110 contains a plurality of source fibers 120 and detector fibers 122 for transmitting 
and receiving light energy, respectively. Each source fiber 120 forms a source detector 
pair with each detector fiber 122 in the imaging head 1 10 to create a plurality of source 
detector pairs. The optical energy entering the target 1 16 at one location is scattered and 
may emerge at any location around the target 116. The emerging optical energy is 
collected by detector fibers 122 mounted in the imaging head 1 10. 

The detector fibers 122 carry the emerging energy to detectors 112. The detectors 
1 12 measure the energy density of the collected optical energy and generate a 
corresponding signal. The detectors may comprise any known photodetector such as 
avalanche photodiodes (APD), silicon PIN photodiodes, silicon photodiodes (SiPD), 
charged coiq)Ie devices (CCD), charge inductive device (CE>), photomuli^lier tubes 
(PMT), multi-channel plate (MCP) and the like. The energy density is, for situations in 
which propagating energy is any type of electromagnetic radiation, equivalent to the 
intensity of said radiation. 
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The data acquisition board 1 14 receives the signal, separates it by wavelength, 
and samples and holds the separated signals for deliver to computer 102. The computer 
102 in turn reads and stores the signal for use in image reconstruction and other analysis. 
Although the above described systems uses DC measurement techniques with an 
5 near infrared energy source, it will be ^predated by those skilled in the art that similar 
systems can be implemented using other measurement techniques, such as time resolved 
measurements and frequency domain methods. In addition measures of acoustic signals 
produced in response to optical absorption (i.e., the photo acoustic effect) as well is use 
of acoustic sources can be similarly implemented. 

10 Data Collection, Data Analysis and Image Reconstruction 

In prior methods for imaging of scattering media, measured data was collected for 
a target medium and an image of the cross-sectional properties of the target was 
generated based on the measured data. The present invention improves on this basic 



0 

\^ process by (1) collecting a time series of measured data and (2) analyzing either the raw 



^ t 15 measured data or a time series of reconstructed images to extract information relating to 
the dynamics of the target medium. 



Data Collection 

As discussed above, previously known systems and techniques sought to collect 
either a time averaged measure th^eby minimizing the effect of dynamic behavior or to 
20 collect a snapshot of the target. Both of these techniques fail to recognize the value of 
information contained in the time domain. 
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The present invention recognizes the value in measuring dynamic properties in 
highly scattering media and provides methods for collecting and processing the measured 
data to generate images of dynamic properties. These methods rely on the high speed 
collection of data using a system such as that described above so that a time series of data 
sets can be collected in the time domain. The time series analysis techniques, discussed 
further below, can then be used to extract the dynamic properties of the target from either 
the measured raw data directly or from a time series of images reconstructed from the 
measured data. 

Data is collected in the present invention using high speed imaging systems such 
as the system illustrated in FIG. 1 . By way of example, referring to FIG. 1, a measured 
data set is obtained by delivering optical energy to each of the source fibers 120 
sequentially and measuring the emerging optical energy at each detector fiber 122 in 
parallel for each sequential source location. This process is repeated over a period of 
time so that a time series of data sets arc generated, each data set represent a complete set 
of data (e.g., detector reading for each source position) for reconstruction of a cross- 
sectional image at an instant in time. 

Reconstruction 

As discussed above, the time series analysis techniques may be employed prior to 
image reconstruction (i.e,, on the raw measured data) or after image reconstruction (i,e., 
on the time series of reconstructed images). Where time series analysis methods are 
applied to raw measured data, image reconstruction will follow based on the processed 
and extracted dynamic information. For example, one method may compute the discrete 
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Fourier transfonn of the measured data and then reconstruct images based on the Fourier 
coefficients at, for instance, a selected frequency. 

Of the known reconstruction methods, one group of methods often adopted to 
analyze the measured data are perturbation methods. Exemplary techniques are disclosed 
in detail in the Barbour '355 patent and the Barbour 4149PC2 application. Briefly, in 
model-based techniques a perturbation formulation based on the radiation transport 
equation, or its approximation the diffusion equation, is generated. The radiation 
transport equation is a mathematical expression describing the propagation of a particle 
through a medium from a source location to a detector location as a function of the 
properties of the medium. The perturbation formulation of this equation relates a 
perturbation in the measure energy density at a detector to a corresponding perturbation 
in the properties of the medium. 

The standard perturbation equation has the following forms: 

W-^x = ^I (1) 

In equation (1), 51 is a vector of source-detector pair intensity differences between 
the measured target medium and the known reference medium (i.e., 51 = I - 1^). W is the 
weight matrix describmg the influence that each volume element CVoxel") of the 
reference medium has on energy traveling from each source to each detector, for all 
source-detector pairs. The volume elements are formed by dividing a slice of the 
reference medium into an imaginary grid of contiguous, non-overlapping pieces. 
Physically, the weight matrix contains the first order partial derivatives of the detector 
responses with respect to the optical coefficients of each volume element of die reference 
medium. 5x is the vector of differences between the known optical properties (e.g., 
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absorption and scattering coefficients) of each volume element of the reference medium 
and the corresponding unknown optical properties of each volume element of the target 
medium. 

A modification to this standard perturbation equation is presented and disclosed in 
the Barbour 4149PC2 ^plication. This modified perturbation equation employs relative 
measurements of energy at the target detectors (e.g., the relative difference between an 
instantaneous measure and a time averaged mean of measures at a detector). This 
method improves the stability of the solution to the perturbation equation and reduces 
sensitivity of the equation to the selected reference medium. 

The modified standard perturbation equation replaces 51 in the standard 
perturbation equation with a proportionate relative difference between two measured 
values multiplied by a reference term of the required units as set forth in the equation (2) 
below: 

^ri^}h (2) 
In equation (2), Ir is the computed detector reading corresponding to a source- 
detector pair of a selected reference medium, and I and lo represent two data 
measurements for a corresponding source-detector pair on one or more targets (e.g., 
background vs. target, or time-averaged mean vs. a specific time point, etc.). The 
resultant term 5Ir therefore rq)resents a proportionate relative difference between two sets 
of measured data. This modified equation limits the effects of modeling errors and 
minimize ill-conditioning of the inverse problem while retaining the correct xmits in the 
solution. 

The corresponding perturbation equation using this modified term is: 
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r ^ r 



(3) 



In equation (3) W, and 6x are the same as W and 5x in equation (1). Referring to 
equations (2) and (3), it can be seen that in the limit, when Ir equals to lo, this equation 
reduces to the standard perturbation formulation shown in equation (1). This formulation 
represents the Bom approximation fomiulation of the modified perturbation equation. A 
similar expression may be written for the Rytov approximation in the following foim: 



The inventive operation performed by equation (2) amounts to a projection 
procedure by which the information content of a measured proportionate relative data 
vector is projected onto a selected reference medium. Apart from simplifying 
measurement requirements, the method represented by equations (3) and (4) also reduces 
susceptibility of the perturbation equation to boundary and optical property variation 
between the target and the reference medium. Additionally, iterative solutions of 
equations (3) and (4) can be easily implemented. In principle, the techniques used in the 
modified perturbation equation can be used to evaluate any measured proportionate 
relative quantity to a reference quantity, for example acoustic signals. 

Time Series Analysis 

As discussed above, time series analysis methods may be applied to either the 
tune series of measurement data or time series of reconstructed images. Applying these 
methods, the present invention enables investigation of the temporal spatial dynamics of 
a target medium. 




(4) 
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Time series analysis techniques are well-known to provide measures of the 
frequency structure and time correlation of time varying processes. In addition, hybrid 
methods such as time frequency analysis (e.g., short time Fourier transform, wavelet 
analysis) can be used to reveal more complex behaviors in physiological systems such as 
the known frequency modulation of the cardiac rate upon respiration. Other linear 
methods can be applied to time domain data such as principle component analysis to 
reveal the time evolution of specific processes in image maps. Practical examples of 
these methods as shown in FIGS, 2-19 can be used to identify, in the case of frequency 
analysis, specific features of the vascular tree, in particular, in the case of imaging human 
tissue, the major arteries in the forearm. This is possible because of the known structure 
dependent frequency response of the vasculature (eg., major arteries beat at a cardiac 
frequency). Further, time correlation methods may be employed to isolate anatomical 
landmarks having known phase relationships (eg., action of antagonistic muscle groups). 
Additionally, similar procedures can be used on the resultant spatial map to produce the 
relevant orthogonal information associated with these functions (i.e., maps of the 
amplitude or phase vs. frequency, correlation vs. time lag, etc.). In addition, complex 
behaviors in the form of nonlinear chaotic activity is often observable in measures of 
physiological systems. As is shown in FIGS. 15-19, the described invention is capable 
of detecting such signatures in a cross-sectional view, in particular as it relates to chaotic 
behavior of the vasculature as measured by time variations in hemoglobin states. While 
such behavior has been observed in surface measurements of tissue, with the exception of 
the instant invention, no other imaging technique is capable of generating equivalent 
information. 
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Experimental Validation 

The following discussion and some of the above examples presents results 
validating the methods and advantages of the present inventions. These examples are 
presented merely as an illustration of the benefits of the present invention and in 
particular of the ability of dynamic optical tomography to provide high contrast images of 
time varying features in dense scattering media. 
Contrast and Resolution Enhancement 

As disclosed above, application of these time-series analysis methods can extract 
significantly greater contrast and resolution from the optical imaging data. Image 
resolution and contrast are important factors for any imaging method. Experience with 
optical tomography has shown that recovered images have low resolution and contrast. 
As an example, refenring to FIG. 2, tomographic data was collected from a laboratory 
vessel measuring approximately 8 cm in diameter and containing four balloons 10-1 
through 1 0-3, two of which are static and the other two were modulated at a fixed 
frequency, with a relative phase of either 0*" or 1 80*. The balloons 10 were suspended in 
a background scattering medium 12 consisting of 2% (v/v) Intralipid and were filled with 
a dilute hemoglobin solution (50 jiM). The balloons 10 were arranged in the vessel 14 as 
shown in FIGS. 2A and 2B. An example of a reconstructed cross-sectional image of 
absorption obtained at a single time point is shown in FIG. 3. Inspection of FIG. 3 shows 
that whereas the four objects structure is detectable, the individual balloons are not well 
resolved, and the recovered contrast is poor. In contrast to this image quality, the map 
shown in FIG. 4, revealing the discrete Fourier transform (DFT) amplitude at the 0.2 Hz 
modulation frequency has nearly 100 fold increased contrast. Also evident is that the 
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spatial resolution of the DFT map is considerably improved compared to the maps shown 
in FIG. 3. The phase portions of the DFT's at 0.2 Hz for the in-phase and the 1 80^-out- 
of-phase cases are shown in FIGS. 5 and 6. The phase relations between the two 
modulated balloons is correctly recovered in each case. 

A similar experiment resolving internal structures based on differences in 
modulation frequency is shown in FIGS. 7 A and 7B. In FIG. 7 A two balloons were 
present, one modulated at 0.1 Hz, the other at 0.24 Hz. Inspection of the reconstructed 
images shown in FIG. 7B, at the corresponding modulation frequencies shows complete 
spatio-temporal resolution. 

Dynamic Physiological Imaging Studies on the Human Forearm 

Following a protocol similar to that used to investigate dynamic behavior in the 
laboratory vessel, the information retrievable from dynamic studies performed on the 
human forearm at rest and in response to specific provocations are also explored. 

Arm at Rest. 

The natural occurrence of vascular frequencies due to respiratory and cardiac 
activity can be exploited to produce a spatial map revealing the presence of different 
functional components of the vascular tree. FIG. 8 shows a map of the logarithm of the 
ratio of the computed DFT amplitudes at the cardiac and respiratory frequencies obtained 
from a times-series of measurements on the forearm at rest. FIG. 9 is a representative 
MR image in the same region of the forearm. An overlay of the two maps having the 
same orientation is shown in FIG. 10. Inspection reveals that in the vicinity of the radial 
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1, interosseous artery 3 and ulnar artery 5, the ratio of the DFT amplitudes (cardiac to 
respiratory) is nearly ten times larger than it is in other regions. The radius 2 and ubia 4 
are also shown. This response can be seen more clearly in FIG. 1 1 A-1 IC, which shows 
the amplitude of the cross-spectral density (CSD) between a surface detector and specific 
5 locations in the image. The particular spectra shown were obtained from points in the 
image corresponding to locations in the flexor digitorum superficialis muscle, and points 
near the radial and interosseous arteries. Inspection clearly reveals that in muscle, which 
contains a dense matrix of microvascular structures, the dominant signal coincides with 
the respiratory frequency. In contrast, we observe a dominant cardiac frequency in the 
10 vicinity of the major arteries. These findings coincide well with the known 
hemodynamic behavior of the different structures in the vascular tree. 

Finger-Flex Study, 

In this study we explored further our ability to measure dynamic behavior in 

1 5 highly scattering media, such as the human forearm, by examining a time series of 

images obtained firom measurements performed on a subject while conducting a finger- 
flex exercise. Finger flexing involves the action of so-called antagonistic muscle groups 
that are located on opposite sides of the forearm, specifically, the flexor digitoium 
superficialis on the ventral side and the extensor digitorum on the dorsal side. Results 

20 shown in FIG. 12 show a map of the amplitude of the DFT at the finger-flexing 

frequency (-0.25 Hz). FIG. 13 shows an overlay of this image onto an MR image of the 
same forearm in the same orientation. Inspection reveals that positions of maximum 
amplitude for fmger-flexing coincide well with the two involved muscle groups. Further 
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evidence supporting the accuracy of this assignment is shown in FIG. 14. Plotted are 
time-series values of recovered absorption values at points in the image coinciding with 
the involved muscles. Noteworthy is the observation that the two signals are 
approximately 180* out of phase with each other, which is the expected response from 
the action of antagonistic muscle groups. 

Imaging the Complexity of the Vascular Response. 

One feature of the disclosed methods that may have considerable practical 
significance is the potential to obtain rehable estimates of the complexity of the vascular 
response. The significance of such measures is that they can reveal features of an 
integrated physiological state that are simply not measurable by using linear methods. 
Accordingly, in this study we have sought to confirm previous reports in the hterature 
that indicate nonlinear chaotic behavior in vascular rhythms, as has been measured 
typically by surface laser Doppler methods. Specifically, we have computed certain 
measures that can reveal nonlinear behavior from time-series image data of the foreann 
of a subject performing a series of deep-breathing exercises. The influence of a 
respiratory stimulus was chosen as a simple noninvasive means to amplify the natural 
oscillatory activity of the vasculature. FIG. 15 shows an example of the magnitude of the 
signal response caused by deep breathing. The measures we have computed from the 
pixel data include k the correlation dimension, and, the largest Lyapunov exponent. 
Positive values for the latter indicate divergence of trajectories in state space, a finding 
consistent with a chaotic time series. A value of zero indicates no net divergence and is 
observed with stochastic data, while negative values indicate system damping. 
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The computed values displayed in a cross-sectional shown in FIG. 16, below 
were obtained using the method of Rosenstein et al (Physica D, 65, 1 17-134 (1993)) for 
small data sets. The image series comprised 240 consecutive time points reconstructed 
from data collected at a rate of 2.8 Hz. The respiration rate was 0.08 Hz for the first 150 
5 time points, and then doubled for the remainder of the measurement period. To assist in 
timing and repeatability, the volunteer was asked to adjust his respiratory rate to match 
the beat of a metronome. 

Referring to FIG. 16, each pixel in the derived image series was band-pass 
filtered over a range that included the vasomotor and respiratory frequencies (0.05-035 

10 Hz). Such filtering colors the noise and can lead to spurious low values for the correlation 
dimension. For this reason data obtained fix>m a stochastic time series in a similar 
fashion was treated similarly. The control data set was obtained by computing an image 
time series fix>m measurements performed on a solid white plastic rod composed of white 
Defrin®, whose scattering prop^es are similar to those of tissue. The rod diameter was 

15 9 cm. 

As shown in FIG. 16, the range of correlation dimrasions computed from the 
physiological image time series is between 2 and 4 for most pixels. Significantly, these 
values are in close agreement with previous studies that measured the correlation 
dimensions of arterial and arteriolar vasomotion in in vitro rat or rabbit preparations. A 
20 similar analysis on the rod data, shown in FIG. 17, reveals significantly higher values, a 
finding consistent with the character of the time series. 

The image in FIG. 18 shows a spatial map of A.\ computed for each pixel. 

Because a stochastic time series can also yield positive values for /li, we have computed 
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surrogate data sets for each pixel and perfonned a test for statistical significance. 
Surrogate data were computed by randomizing the phase of the Fourier spectrum 
followed by recovery of the resultant time series from which the surrogate A\ were 
computed. FIG. 19 shows a map of the significance level, in each pixel, for the 
5 physiological data compared to surrogate data. Inspection reveals that most areas of the 
map have positive Jii values that fail the null hypothesis. Note that the optical measures 
performed are known to exhibit considerable sensitivity to the microvessels. These 
findings, indicating the occurrence of nonlinear chaotic activity of the vasculature 
revealed in a cross-sectional view, thus coincide well with previous reports whose 
1 0 measures were limited to the study of vascular dynamics in the near surface of tissue. Its 
worth emphasizing that the above findings are the first of their kind and, to the best of our 
knowledge, no other imaging modality is capable of providing equivalent measures of 
vascular dynamics. 

Strategies of the types explored can be used to investigate a range of vascular 
15 anomalies. For example, the deep breathing exercise may aid in the diagnosis of cancer. 
Solid tumors have altered vascular beds. Enhanced breathing will cause redistribution of 
the flow patterns in tissue. Regions having altered vascular architecture should behave 
very different firom those with a normal architecture. A similar approach can be used to 
identify tissue regions "at risk" due to the effects of diabetes, atheroscleorsis (small 
20 vessel disease), or long term effects of smoking (large vessel disease). 

Although the examples above focus on near infrared energy sources for imaging 
human tissue, the methodology of the present invention is applicable to essentially any 
wavelength for any energy source (e.g., electromagnetic, acoustic, etc.), responsive to 
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dynamic properties of a scattering medium and for any source condition (e.g., time 
independent, time harmonic, time resolved). Moreover, although the examples above 
investigate cross-sectional imaging, it will be appreciated that these methods are 
applicable to and include volumetric measures and imaging. 

Further, although illustrative embodiments have been described herein in detail, 
those skilled in the art will appreciate that variations may be made without departing 
from the spirit and scope of this invention. Moreover, unless otherwise specifically 
stated, the terms and expressions used herein are terms of description and not terms of 
limitation, and are not intended to exclude any equivalents of the system and methods set 
forth in the following claims. 
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1 . A method for enhanced imaging of a target medium comprising: 
directing energy into a target medium from at least one source during a period of 

time, the target medium having dynamic properties during the period of time; 

measuring the density of the energy emerging from the target medium during the 
period of thne using at least one detector, the density of the energy emergmg from the 
target medium being a ftmction of the dynamic properties of the target medium; and 

generating a map of the dynamic properties of the target medium based on the 
measured density of energy emerging from the target medium. 

2. The method of claim 1, ftuther comprising: 

generating a time series of images of the properties of a target medium based on 
the measured energy emerging from the target medium, wherein each image represents 
the cross-sectional properties of the target medium at a time interval during the period of 
time, 

3. The method of claim 2, wherein generating the map of the dynamic 
properties of the target medium comprises processing the time series of images using 
time s^es analysis methods. 

4. The method of claim 1 , wherein generating the map of the dynamic 
properties of the target medium comprises processing the measured energy at each 
detector using time series analysis methods. 
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5. The method of claim 1, wherein the image of the dynamic properties is 
generated using time series analysis methods. 

6. The method of claim 5, wherein the time series analysis methods comprise 
linear time series methods. 

?• The method of claim 6, wherein the linear time series analysis methods are 
at least one of frequency analysis, time conrelation analysis, time frequency analysis and 
principle component analysis. 

8. The method of claim 5, wherein the time series analysis methods comprise 
nonlinear time series analysis methods. 

9. The method of claim 1 , wherein measuring the energy comprises a 
plurality of measurements collected at a sampling rate not less than twice the reciprocal 
of the highest frequency of a dynamic property to be imaged. 

10. The method of claim 1, fiirther comprising applying a provocation to the 
target medium. 

1 1 . The method of claim 1 , wherein the target medium is human tissue having 
a vascular tree and the provocation has a dynamic effects on the vascular tree. 

12. The method of claim 11, wherein the provocation is an autonomic 
stimulus. 

1 3. The method of claim 1 1 , wherein the provocation is a local stimulus. 

14. The method of claim 1, wherem the energy is optical energy having a 

wavelength in the near infrared region of the electromagnetic spectrum. 
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15. The method of claim 14, wherein the optical energy directed toward the 
medimn includes at least two wavelengths of near infrared energy. 

16. The method of claim 14, wherein the target medium is human tissue 
having a vascular tree containing blood, the vascular tree comprising veins, arteries and 
micro vessels, the blood having time varying absorption and scattering properties to the 
near infrared energy as a function of blood oxygenation and blood volume. 

17. The method of claim 16, wherein generating an image of the dynamic 
properties of the target medium comprises generating an image of at least one of the time 
varying absorption and scattering properties of the target medium. 

1 8 . The method of claim 1 6, wherein generating an image of the dynamic 
properties of the target medium further comprises using time series analysis to enhance 
the contrast of at least one of veins, arteries and micro vessels. 

1 9. The method of claim 1 , wherein the energy is optical energy in the visible 
spectrum. 

20. The method of claim 1 , wherein generating an image of the dynamic 
properties of the target medium further comprises processing the measured energy using 
a modified perturbation formulation of the radiation transport equation, wherein the 
modified perturbation formulation uses relative energy measurements. 

2 1 . The method of claim 20, wherein the relative energy measurements are the 
relative differences between a measure at an instant in time and a time average mean of 
measxires during the period of time, 
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22. A system for enhanced imaging of a target medium comprising: 

a source for directing energy into a target medium from at least one source during 
a period of time, the target medium having dynamic properties during the period of time; 

a detector for measuring the density of the energy emerging from the target 
medium during the period of time using at least one detector, the density of the energy 
emerging from the target medium being a function of the dynamic properties of the target 
medium; 

a data acquisition means for receiving the measured energy densities during the 
period of time; and 

a computer connected to the data acquisition means, the computer having code for 
generating a map of the dynamic properties of the target medium based on the measured 

density of energy emerging from the target medium. 

23. The system of claim 22, v/herein the computer further includes code for 
generating a time series of images of the properties of a target medium based on the 
measured energy emerging from the target medium, wherein each image represents the 
cross-sectional properties of the target medium at a time interval during the period of 
time. 

24. The system of claim 23, wherein the computer further includes code for 
processing the time series of images using time series analysis methods to generate the 
map of the dynamic properties. 
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25. The system of claim 24, wherein the computer fiirther includes code for 
processing the measured energy at each detector using time series analysis methods to 
generate the m^ of the dynamic properties of the medium. 
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Docket No. 0887-4 147PC2 



COMBINED DECLARATION AND POWER OF ATTORNEY FOR 
ORIGINAL, DESIGN, NATIONAL STAGE OF PCT, SUPPLEMENTAL 
DIVISIONAL, CONTINUATION OR CONTINUATION-IN-PART APPLICATION 

As a below name inventor, I hereby declare that: 

My residence, post office address and citizenship are as stated below next to my name, 

I believe I am the original, first and sole inventor (if only one name is listed below) or an original, first and joint 
inventor (if plural names are listed below) of the subject matter which is claimed and for which a patent is sought on 
the invention entitled: 

METHOD AND SYSTEM FOR IMAGING THE DYNAMICS OF A SCATTERING MEDIUM 

the specification of which 

a. [ ] is attached hereto 

b. [ ] was filed on as application Serial No. and was amended on 

. (if appUcable). 

PCT FILED APPLICATION ENTERING NATIONAL STAGE 

c. [x] was described and claimed in International Application No. PCT/USQQ/25136 filed on 14 

September 2000 and as amended on . (if any). 

I hereby state that I have reviewed and understand the contents of the above-identified specification, including the 
claims, as amended by any amendment referred to above. 

I acknowledge the duty to disclose information which is material to the patentability as defined in Title 37, Code of 
Federal Regulations, § 1 .56. 

I hereby specify the following as the correspondence address to which all communications about this application are 
to be directed: 

SEND CORRESPONDENCE TO: MO RGAN & FINNEGAN, L.L.P 

345 Park Av enue 
New York, N .Y. 10154 

DIRECT TELEPHONE CALLS TO: Tod M. Melgar 
(212)758-4800 

[ j I hereby claim foreign priority benefits under Title 35, United States Code § 1 19(a)-(d) or under 
§ 365(b) of any foreign application(s) for patent or inventor's certificate or under § 365(a) of any PCT international 
application(s) designating at least one country other than the U.S. listed below and also have identified below such 
foreign application(s) for patent or inventor's certificate or such PCT international application(s) filed by me on the 
same subject matter having a filing date within twelve (12) months before that of the application on which priority is 
claimed: 

[ ] The attached 35U.S.C. § 119 claim for priority for the application(s) listed below forms a part of this 
declaration. 
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Application Date of filing Date of Issue Priority 

Country/PCT Number (day, month, vr^ (day, month, vr^ Claimed 

r 1 YES r 1 NO 

r 1 YES r 1 NO 

[ ] YES r 1 NO 



[X] I hereby claim the benefit under 35 U.S.C. § 119(e) of aay U.S. provisional application(s) listed below. 
Provisional Application No. Date of Filing Tdav. month, vr^ 

60/153,926 14 September 1999 

60/154,099 15 September 1999 



ADDITIONAL STATEMENTS FOR DIVISIONAL, CONTINUATION OR CONTINUATION-IN-PART 
OR PCT INTERNATIONAL APPLICATIONfSl (DESIGNATING THE U.S.^ 

I hereby claim the benefit under Title 35, United States Code § 120 of any United States application(s) or under 
§ 365(c) of any PCT intemational application(s) designating the U.S. listed below. 

US/PCT Application Serial No. Filing Date Status (patented, pending, abandoned)/ 

U.S. application no, assigned (For PCT) 



US/PCT Application Serial No. Filing Date Status (patented, pending, abandoned)/ 

U.S. application no. assigned (For PCT) 

[ ] In this continuation-in-part application, insofar as the subject matter of any of the claims of this 
application is not disclosed in the above listed prior United States or PCT intemational appiication(s) in the manner 
provided by the fu-st paragraph of Title 35, United States Code, § 1 12, 1 acknowledge the duty to disclose material 
information as defined in Title 37, Code of Federal Regulations, § 1 .56(a) which occurred between the filing date of 
the prior application(s) and the national or PCT intemational filing date of this application. 

I hereby declare that all statements made herein of my own knowledge are true and that all statements made on 
information and belief are believed to be true; and further that these statements were made with the knowledge that 
willful false statements and the like so made are punishable by fme or Imprisonment, or both, under Section 1001 of 
Title 18 of the United States Code and that such willful false statements may jeopardize the validity of the 
application or any patent issued thereon. 
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I hereby appoint the following attorneys and/or agents with full power of substitution and revocation, to prosecute 
this application, to receive the patent, and to transact all business in the Patent and Trademark Office' connected 
therewith: John A. Diaz (Reg. No. 19,55.0), John C. Vassil (Reg, No. 19,0981 Alfred P. Ewert (Reg No UM2) 
David H. Pfeffer (Reg. No. n^SZl), Harry C. Marcus (Reg. No. 223fi), Robert E. Paulson (Reg. No. 21,046 ) 
Stephen R. Smith (Reg. No. 22,615), Kurt E. Richter (Reg. No. 24,0521 J. Robert Dailey (Reg. No. 27,434 ) Eugene 
Moroz (Reg. No. 2X222), John F. Sweeney (Reg. No. 27.471\ Arnold 1. Rady (Reg. No. 26.601 1 Christopher A 
Hughes (Reg. No. 26,914), William S. Feiler (Reg. No. 26.7281 Joseph A. Calvaruso (Reg. No, ZSOSJ) James W 
Gould (Reg. No. 28^), Richard C. Komson (Reg, No. 27,913 1 Israel Blum (Reg. No. l^JUiX Bartholomew 
Verdirame (Reg. No, 28,483), Maria C.H. Lin (reg. No. 22*223), Joseph A. DeGirolamo (Reg. No. 28.5951 Michael 
P. Dougherty (Reg. No. 32J2£Ll Seth J. Atlas (Reg. No. 32.454 ). Andrew M. Riddles (Reg. No. 3C657 ) Bruce D 
DeRenzi (Reg. No.,33,676}, Michael M. Murray (Reg. No. 32,5371 Mark J. Abate (Reg. No. 32,5271 John T 
Gallagher (Reg. No. U^), Steven F. Meyer (Reg. No. 15,613), Kenneth H. Sonnenfeld (Reg. No.Ji2S5) Tony 
V. Pezzano (Reg, No. I8,22igu Andrea L. Wayda (Reg. No . 43.979 ;^ and Walter G. Hanchuk Reg. No. (35,179 ) of 
Morgan & Finnegan, L.L.P. whose address is: 345 Park Avenue, New York, New York, 10154; and Michael S 
Marcus (Reg. No. 31,727) and John E. Hoel (Reg. No. 26.279^ of Morgan & Finnegan, L.L.P., whose address is 
1775 Eye Street, Suite 400, Washington, D.C. 20006. 

[ ] I hereby authorize the U.S. attorneys and/or agents named hereinabove to accept and follow instructions 
from 



as to any action to be taken in the U.S. Patent and Trademark Office 



regarding this application without direct communication between the U.S. attorneys and/or agents and me. 
In the event of a change in the person(s) from whom instructions may be taken I will so notify the U.S. 
attorneys and/or agents hereinabove. 

Full name of sole or iHrst inventor Randall L. BARBOUR 

Inventor' s signature * 

^ date' 
Residence 1 5 Cherry Lane. Glen Head, New York 1 1 546 My 




Citizenship United States 



Post Office Address same as above 



Full name of second joint inventor, if any 

Inventor's signature* 

date 

Residence 

Citizenship 

Post Office Address 
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Full name of third joint inventor, if any 

Inventor's signature* 

date 

Residence 

Citizenship 

Post Office Address 



[ ] ATTACHED IS/ARE ADDED PAGE(S) TO COMBINED DECLARATION AND POWER OF 
ATTORNEY FORM FOR SIGNATURE BY FOURTH AND SUBSEQUENT INVENTORS 

* Before signing this declaration, each person signing must: 

1 . Review the declaration and verify the correctness of all information therein; and 

2, Review the specification and the claims, including any amendments made to the claims. 

After the declaration is signed, the specification and claims are not to be altered. 

To the inventor(s): 

The following are cited in or pertinent to the declaration attached to the accompanying application: 

Title 37. Code of Federal Regulation, § 1.56 

Duty to disclose information material to patentability. 

(a) A patent by its very nature is affect with a public interest. The public interest is best served, and 
the most effective patent examination occurs when, at the time an application is being examined, the Office 
is aware of and evaluates the teachings of all information material to patentability. Each individual 
associated with the filing and prosecution of a patent application has a duty of candor and good faith in 
dealing with the Office, which includes a duty to disclose to the Office all information known to that 
individual to be material to patentability as defined in this section. The duty to disclose information exists 
with respect to each pending claim until the claim is canceled or withdrawn jfrom consideration, or the 
application becomes abandoned. Information material to tne patentability of a claim that is canceled or 
withdrawn from consideration need not be submitted if the information is not material to the patentability 
of any claim remaining under consideration in the application. There is no duty to submit information 
which is not material to the patentability of any existing claim. The duty to disclose all information known 
to be material to patentability is deemed to be satisfied if all information known to be material to 
patentability of any claim issued in patent was cited by the Office or submitted to the Office in the manner 
prescribed by §§ 1 .97(b)-(d) and 1 .98. However, no patent will be granted on an application in connection 
with which fi'aud on the Office was practiced or attempted or the duty of disclosure was violated through 
bad faith or intentional misconduct. The Office encourages applicants to carefully examine: 

(1) prior art cited in search reports of a foreign patent office in a counterpart application, and 
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(2) the closest information over which individuals associated M^ith the filing or prosecution of 
a patent application believe any pending claim patentably defines, to make sure that any 
material information contained therein is disclosed to the Office. 

Title 35. U.S. Code § 101 
Inventions patentable 

Whoever invents or discovers any new and useful process, machine, manufacture, or composition of 
matter, or any new and useful improvement thereof, may obtain a patent therefor, subject to the conditions and 
requirements of this title. 



Title 35 U.S. Code ^ 102 

Conditions for patentability; novelty and loss of right to patent 
A person shall be entitled to a patent unless - 

(a) the invention was knovm or used by others in this country, or patented or described in a printed 
publication in this or a foreign country, before the invention thereof by the applicant for patent, 

(b) the invention was patented or described in a prmted pubHcation in this or foreign country or in 
public use or on sale m this country, more than one year prior to the date of application for patent in the United 
States, or 

(c) he has abandoned the invention, or 

(d) the invention was first patented or caused to be patented, or was the subject of an inventor' s 
certificate, by the applicant or his legal representatives or assigns in a foreign country prior to the date of the 
application for patent in this country on an application for patent or inventor's certificate field more than twelve 
months before the filing of the application in the United States, or 

(e) the invention was described in a patent granted on an application for patent by another filed in the 
United States before the invention thereof by the applicant for patent, or on an international application by another 
who has fulfilled the requirements of paragraphs (1), (2), and (4) of section 371(c) of this title before the invention 
thereof by the applicant for patent, or 

(f) he did not himself invent the subject matter sought to be patented, or 

(g) before the apphcant's invention thereof the invention was made in this country by another had not 
abandoned, suppressed, or concealed it. In determining priority of invention there shall be considered not only the 
respective dates of conception and reduction to practice of the invention, but also the reasonable diligence of one 
who was first to conceive and last to reduce to practice, fi-om a time prior to conception by the other 

Title 35. U.S. Code 6 103 

Conditions for patentability; non-obvious subject matter 

A patent may not be obtained though the invention is not identically disclosed or described as set forth in 
section 102 of this title, if the differences between the subject matter sought to be patented and the prior art are such 
that the subject matter as a whole would have been obvious at the time the invention was made to a person having 
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ordinary skill in the art to which said matter pertains. Patentability shall not be negatived by the manner in which 
the invention was made. 

Subject matter developed by another person, which qualifies as prior art only under subsection (f) or (g) of 
section 102 of this title, shall not preclude patentability under this section where the subject matter and the claimed 
invention were, at the time the invention was made, owned by the same person or subject to an obligation of 
assignment to the same person. 



Title 35. U.S. Code § 1 12 Tin p art^ 



Specification 



The specification shall contain a written description of the invention, and of the manner and process of 
making and using it, in such full, clear, concise and exact terms also enable any person skilled in the art to which it 
pertains, or with which it is mostly nearly connected, to make and use the same, and shall set forth the best mode 
contemplated by the inventor of carrying out his invention. 

Title 35. U.S. Code 6 119 

& 



|:j Benefit of earlier filing date in foreign country; right of priority 



m 



An application for patent for an invention filed in this country by any person who has, or whose legal 
^„ representatives or assigns have, previously regularly filed an application for a patent for the same invention in a 

foreign country which affords similar privileges in the case of applications filed in the United States or to citizens of 
^ the United States, shall have the same effect as the same application would have if filed in this country on the date 
? on which the application for patent for the same invention was first filed in such foreign country, if the application in 
Cl this country is filed within twelve months fi-om the earhest date on which such foreign application was filed; but no 
y patent shall be granted on any application for patent for an invention which had been patented or described in a 
yk. prmted publication in any country more than one year before the date of he actual filing of the application in this 
^'j country, or which had been in public use or on sale in this country more than one year prior to such filing. 



nj Title 35. U.S. Code S 120 



Benefit or earlier fihng date in the United States 



An application for patent for an invention disclosed in the manner provided by the first paragraph of section 
1 12 of this title in an application previously filed in the United States, or as provided by section 363 of this title, 
which is filed by an inventor or inventors named in the previously filed application shall have the same effect, as to 
such invention, as though filed on the date of the prior application, if filed before the patenting or abandonment of or 
termmation of proceedings on the first application or an application similarly entitled to the benefit of the filing date 
of the first application and if it contains or is amended to contain a specific reference to the earlier filed appHcation. 

Please read carefully before signing the Declaration attached to the accompanying AppUcation. 

If you have any questions, please contact Morgan & Finnegan, L.L.P. 



FORMrCOMB-DEC.NY 
Rev. 10/00 
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